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A Dynamic Model for 
Information Incorporating the Effects 
of Counter-CJ 



Paul H. Moose 



I . Introduction 

In modern warfare, an operational commander is intimately concerned 
with the quality, timeliness and completeness of his "picture" of the 
tactical situation. To a very large extent, his fortunes and those of 
his assigned forces depend on his having available, when and where he 
needs it, accurate data about the status, location and activities of both 
his own and the enemies' forces. Similar requirements extend well down 
into subordinate echelons of his command, including individual unit 
commanders and even "smart" weapons. 

The methods by which necessary information is acquired are remarkably 
diverse. Included, in our "computer age", are sophisticated radar and 
intercept equipments, a variety of imaging systems and acoustic sensors 
as well as ordinary direct reports and observations from the commander's 
own personnel. Inputs from special intelligence channels, and many other 
categories of reports constantly arriving by a variety of means round out 
a massive and continuing informational input. The staff, assisted undoubt- 
edly by modern automatic data processing equipments, is regularly creating 
and updating their assessment of the situation in order to give the best 
operational picture they can to their commander. The commander will, to 
a very great degree, make rational and reasonably predictable decisions 
for the future activities of his forces based on the world view he has 
developed from this sequence of images. 



II. A View of Information 



We wish at this point to make several important observations about 
this "image of reality" that the commander works with. First, the images 
he has are never absolutely correct, that is, they contain errors. Nor 
are they perfectly sharp, that is, there are always many questions that are 
unanswered, or elements of contradiction or ambiguity. Secondly, an image 
gets fuzzier and fuzzier the further into the future one attempts to extra- 
polate it. This is because most elements of the picture are dynamic, i.e. 
they change (location, behavior, etc.) with time. The attributes of the 
elements may be partially constrained. For example, a ship cannot move 
faster than about 30 knots. However, after sufficient time most features 
of the picture will have totally relaxed, and may have taken on any of 
their possible values or conditions. 

This "fuzzyness in the crystal ban" axiom has a corollary. If the 
commander loses, or turns off, his senses or sources of information, his 
"current image" will grow fuzzier and fuzzier with time until it is 
completely blurred. Put another way, a commander only maintains his 
uncertainty about what is going on below its worst possible level by 
virtue of the continual application of systemic resources to guarantee an 
inflow of new information. 

The second law of thermodynamics holds that entropy(disorder) will 
grow to its maximum possible value within the constraints of the system. 

We shall identify uncertainty with entropy. Thus sensory devices and 
information sources provides the constraints on uncertainty. They do 
this by continually importing information (negentropy) to offset 
uncertainty's inevitable growth. 
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Ill . Birth and Death of Uncertainty 



We postulate that if left unconstrained, uncertainty will grow from 

its current value toward its maximum possible value H with a rate 

m3 x 

proportional to the remaining knowledge. (H measures the worst case, 

max 

total ignorance, where all possibilities are equally probable.)^ If T s 
is the systemic relaxation time, then 



"b ■ ^ 



(h max- h > 



( 1 ) 



models the birth of uncertainty, Hg, as proportional to the remaining 
knowledge, -H, divided by the system relaxation time T . Left 
unconstrained, uncertainty will grow as shown in Figure 1. 




^Growth of this type is not an unreasonable assumption as can be 
seen in the work of Moose and Harrison, "An Analytic Model of Coordinated 
Effort with Application to Surveillance C^", May 1979, AD # A071 -081 . 



3 



However, the commander is continually receiving new data that helps 
to reduce uncertainty; he is receiving "negentropy" . We assume that the 
rate of entropy reduction is directly proportional to his state of 
ignorance. Thus 

H d = -Tj 1 H (2) 

expresses the fact that arriving information destroys uncertainty at a 
rate in direct proportion to its current value. Tj is a characteristic 
negentropy arrival time. Let us further assume that the causes (1) and 
(2) may be superimposed such that 

S ' T s' "max -< t s' + t ;’> h < 3 > 

models the net growth (or reduction) in uncertainty at any point in time. 
It is clear that this system is stationary when 



H MAX T S + T I 

since at this point the birth and death rates exactly cancel one another 

O 

and H = 0. 

Although this simple first order linear system is not very sophisti- 
cated, we see that it already illustrates some practical features. When 
T<., the systemic relaxation time is very short compared to the character- 
istic arrival time of new information, uncertainty finds its equilibrium 
near maximum. On the other hand, if the relaxation time is equal to the 
characteristic arrival times, uncertainty will be cut to one-half its 
maximum value. When information arrives at a much greater rate than the 
attributes of the picture can change, uncertainty is very low and finds 
equilibrium at a level in direct proportion to the ratio "^I/T^. One must 
be cautious and take due note of the fact that Tj & T $ (or Vj = Tj and 
= TjJ ) are global or macroscopic system variables. Vj & might have 
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units, for example, of bits/bit per hour and refer to the average behavior 
of the entire system ensemble much as species birth and death rates are 
typically measured in births (or deaths)/unit of population per year, and 
describe an average of the entire population. 



IV. Information War 
( 2 ) 

Rona' has described the concept of "information war" as a dominant 

factor in the conduct of modern warfare. In an information war, one 

actively attempts to deny the enemy knowledge of his force positions, 

numbers, intentions, etc. This is done by a variety of means. Included, 

for example, are cover and deception tactics, distribution of radar chaff, 

decoys, false messages, etc. One also works to keep his own communications 

intact and secure, but intercepts and exploits and/or jams those of the 

3 

enemy. One may also try to physically disable enemy C facilities and 

channels. In all of this, the purpose is to try to reduce one's own 

uncertainty by assuring a steady, reliable inflow of relevant information, 

a term we have already described above. But moreover, to disrupt the 

opposition's flow of information and ultimately blur or distort his image 

of the operational situation. This will be cause, we maintain, for poor 

decisions on his part thereby enhancing one's own force effectiveness. 

Let us suppose X & Y represent the entropies of two opposing sides. 
° 3 

We should include, in X, a "counter-C term to represent informational 



disruption by Y, and vice versa. Thus, the rate equations take the form 

(4) 



O 

X 



V SX X MAX " ^ V SX + V IX )X + r YX^ Y MAX" Y) 



Y " V $Y Y MAX " ^ V SY + V IY )Y + r XY^ X MAX " X ^ 



( 2 ) 

Rona , T.P., "Weapon Systems and Information War", Boeing 
Aerospace Co., Seattle, WA, July 1976. 
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We have assumed that Y's disruption of X 1 s flow of information is 

directly proportional to his current knowledge, anc * likewise for 

3 

X's disruption of Y. The "Counter-C coefficients ", r yx & r XY are, 
we presume, positive if in fact the information war is having the desired 
results, at least on the average. 

There is a question, however, about how these coefficients are to be 
chosen. In particular does r YX depend only on Y's efforts against X 
or is it also dependent on X's own knowledge, and if so, how? It seems 
plausible to suppose that if X is very short of knowledge of the situation 
already, he may be difficult to confuse even more, whereas if his knowledge 
is great he may be much more susceptible to disruption, deception, decoys 
and jamming, Admittedly, this is a highly speculative argument, but is is 
an extremely important point because it determines whether the two systems 
are linearly or non-linearly coupled, which, as we shall see presently, has 
an immense influence on their dynamic behavior. 

* c 

Let us list some options for Y coupling to X. 





r YX 




Form 


Result 


a) 


r YX = 


y YX X MAX ^ 0 


constant 


linear coupling 


b) 


r YX = 


y YX^ X MAX " * 0 


proportional to 
X's knowledge 


2nd degree 
(non-1 inear) 










coupl ing 


c) 


-< 

X 

II 


-< 

X 

X 

\V 

o 


proportional to 


2nd degree 




X' s uncertainty 


(non-linear) 
coupl ing 






d) 


r YX = 


y yx x max e ' X * 0 


increases exponen- 
tially with x' s 
knowledge 


non-linear 
coupl ing 


e) 


r YX = 


Y YX x < x max- x > 40 


maximum coupling 


3rd degree 




whe„ „ . X Hflx and 


(non-linear) 








2 


coupl ing 



no coupling when 
X = 0 or X = X HAX 



6 



If we chose case b) above as the most intuitively appealing dependence, 
then Eq's (4) take the form 

X = ^ V SX + y YX Y MAX^ X MAX ‘ ^ V SX + V IX + y YX Y MAX^ X ' y YX Y ^ X MAX ' X ^ 

(5) 

Y = ^ V SY + Y XY X MAX^ Y MAX ' ^ V SY + V IY + Y XY X MAX^ Y " Y XY X ^ Y MAX " Y ^ 
which are 1st order, 2nd degree - non-linear coupled differential equations. 
We may write them more compactly as, 

X = a o X MAX " a l X ' a 2 Y ^ X MAX ' X ^ 

( 6 ) 

Y= g o Y MAX- ^ Y MAX- Y ) 



To summarize, we have the following relationships between coefficients and 
definitions of system parameters: 

A. Coefficient Relationships 



i) 


a o = 


(V S X + a 2 Y l 


MAX^ ’ B o “ ^ V SY + B 2 X MAX^ 


ii) 


a l = 


a o + V I X 


* B 1 = B o + V IY 


i 1 i ) 


a 2 = 


y yx 


b 2 = y xy 


Parameter Definitions 




i) 


X MAX 


* y max ; 


Maximum Uncertainty (Bits) 


ii) 


< 

Ln 

X 


V SY ; 


Uncertainty Birth Rates (bits/bit per unit time) 


i i i ) 


V IX’ 


V IY > 


Uncertainty Death Rates resulting from data 
inputs (bits/bit per unit time) 


iv) 


X 

>- 

>~ 


y xy ; 


3 -1 

Counter C coefficients (bits per unit time) 



O 

Note that if the counter C° coefficients are both zero, Eq's (6) revert to 
a pair of uncoupled equations of the form given by Eq. (3). 
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V. Analysis 



Analysis of the behavior of the non-linearly coupled 1st order 
differential equations, presented in Eq's (6) as a model for information 
dynamics between opposing C systems, is usefully subdivided into separate 
treatments of 1.) "Stationary" or "Equilibrium" conditions, 2.) dynamic 
behavior near equilibrium points, and 3.) dynamic behavior far from 
equilibrium points. Before proceeding to the discussion of each of these, 
it is convenient to consider a normalized version of the equations as 
follows. Let 



x * f (v sx + V ' (v sx + v ix + V X ' Y yx y(, - x)] Y max 
y * ttV * V -'V * 1y + V y - Y xy x(, - y>) X max 



(7) 



where we now have defined normal ized entropy variables and entropy rates, 



X = X/X MAX’ X ~ X/X MAX ’ y = Y/Y MAX’ and y = Y/Y MAX* 

Note that 0 £x, y^ 1 for physically realizable conditions. We may 
again introduce the more compact form 



x = (“o - a,x - a 2 y(l-x) Y max 

y * (6 o ' B l y - 6 2 x(, - y » X max 
with the coefficients defined now as follows: 



( 8 ) 



A. Normalized Coefficients 



0 % * v sx * Y yx 
ii) a, - o + v, 



3 = v + Y 

o sy xy 



3, = 3 + v. 

P 1 p o l y 



iii) a 2 = y 



yx 



Y 



xy 
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B. 



Normalized Parameters 



^ v sx " V SX /Y MAX 



v $y = ^sy^MAX’ Normalized Birth Rates 
(bits'Vunit time) 



ii) ^IX^MAX 9 V IY ^IY^^MAX* Normalized Death Rates 




3 -1 

iii) y , y ; Counter C coefficients (bits - /unit time) 
^ the same as in Eq 1 s (6)). 



V.l Stationary Points & Sensitivity 



Stationary points will occur when x = y = 0, that is when 



the conditions 



a 



o 



a-| x - a£ y(l-x) = 0 



( 9 ) 



6 0 8-j y - S 2 x(l-y) = 0 



are met by x and y. 

Eq (9) may be solved for the values of x and y that yield these 
conditions. Because of the presence of the product term xy, there are, 
in general, two stationary points, which we shall designate (x-j , y^ ) and 
(* 2 * ^2^' before we at the general solutions of (9), consider the 
case where y x ^ = ^ = 0 and y yx = a.^ > 0* This case, which we shall 

designate as "Case 1", models the conditions where X is exercising no 

2 3 

counter-C on Y, but Y is actively counter-C ing X. We wish to see 

what shifts occur in the equilibrium conditions (x Q , y Q ) as a function 

3 

of Y's counter-C effort y 

'yx 

Observe that the equations are no longer quadratic and there is only 
one equilibrium point (x-| , y-| ) and since 3 2 = 0, y - 1 = y Q so that 
immediately we find 



a 



a 



o 
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(10) 



1 s 1 s 

Thus we see that Y state of knowledge is left unchanged but X 

state deteriorates from x Q to x-j , (x-j £ x Q ), as a function of Y 

countering effort = y^ x in accordance with Eq (10). Figure 2 shows 

a plot of X s relative increase in entropy, (X. - x J/x , versus y /v , 

i u o y x s x 

• c 1 c 

Y countering efforts in relation to X natural entropic birth rate. 

Two curves are shown; the lower is for the case x Q = 0.1, y Q = 2/3 and 

the upper is for the case x Q = 0.1, y Q = 0.1. Note they both approach the 

value 9 for large values of y . This is maximum entropy, x, = 1 , 

yx i 

• c 

since we started at x q = .1; i.e. the most Y can do is increase X 

entropy by 9 times since that represents total uncertainty for X. Note 

for small values in Figure 2 that the curves approach straight 45° lines on 

log-log scales indicating that (x-,-x )/x and y, /v / v CY are linearly 

related for small values y . Indeed careful expansion of Eq (10), retaining 

yx 

all first order terms in a 9 , (a ? = y ), shows that for small y . 

l. l. y a y x • 

< X 1 - x o )/x o" ^ • O-X )d-y 0 ). (11) 

sx 0 

Readers may easily verify for themselves that Eq (11) gives the same 

results as are shown in Figure 2 for values of y /v ^1. 

3 'yx sx 

Eq (11) gives the "counter-C sensitivity" of the system for one way 
coupling. It shows that if either X or Y are in an initial state of great 
uncertainty, counter-C efforts of one against the other will be of little 
value since either (l-x Q ) or (1-y ), or both, are approximately zero. 
However, if both sides are functioning with high informational efficiency 
each side is maximally vulnerable to countering efforts by the other. 
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Having obtained an initial feeling for the sensitivity of the 
equilibrium point to the coupling factors, let us return to the general 
problem posed by Eq's (9) where both sides are actively engaged in 
counter-C efforts. Solving for the roots of Eq's (9), we find 

y l ,2 = C 1 X 1 ,2 + c o 

X 1 2 = P^ + Y (P^ " 4q) 2 



( 12 ) 



where 



Ct, 



P = 1 + 



a 2 c l 



(13) 



q = 



a 



a 2 C l 




(14) 



Cq ~ (c^Sq " a 0 ^2^ a 2^1 " ) 


(15) 


C-| — ( ot-j 3 ^ " a 2^2^^ a 2^l" ^2^ 


(16) 



n 

It can be shown that p - 4q > 0 and therefore the roots (x-j , y-j ) and 
(x 2 , y 2 ) ace always real. Further general analysis of Eq's (12) has 
shown that one of the roots is always found in the physically realizable 
space (0 < x-j < 1 , 0 < y^ < 1) and the other is found outside the space, 
but in the first quadrant. ^ 

Considerable practical insight can be obtained by considering the 
following "matched case" which we shall designate "Case 2". For Case 2, 
we let 



lyx . ^2 ^sx 
y xy v sy 



( 17 ) 



( 3 ) 



K.E. Woehler, 



"Root Locations", private communication, 1980. 
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that is we keep the counter-C J efforts in a constant ratio equal to the 
ratio of natural entropic birth rates. With this constraint, we find that 



= 0 , c, = yjx r 



( 18 ) 



'0 ~ ’ ~1 ■' o ' "0 

so that 

y l,2 = (y o /x o^ X 1 ,2 ’ (19 ^ 

For Case 2, we see that both the equilibrium points must always lie on a 
straight line in the x,y plane passing thru the origin and the uncoupled 



stable point x o ,y Q . Also, we note that 



x. a-j 



p ^>Oandq - 'f £ >0 



( 20 ) 



2 ^o "2 

and therefore all roots are positive. In Figure 3 we have plotted the 
relative locations of the roots, (x-j ^ ~ x 0 )/x Q > versus Y yx^ v sx' t ^ 1e 

3 

normalized counter-C coefficient for the initial conditions x = y = 0.1. 

o o 

Note that one root pair starts at (°°, ») and moves towards (1,1) as 

Yy X increases whereas the other root pair begins at (x Q , y Q ) and moves 

toward (1,1) with increasing y (and increasing y according to 

yx xy 

Eq (17)). 

2 ^ 

By a very careful expansion of the radical (p - 4q) 2 with the 
additional constraint that c^=l, i.e. that x Q =y o (met by the initial 
conditions used in Fig. 3), we find that for small y ; 

X 1 " x o Y VX 

“4“* * ft (, - x ) (1 - y ) (21) 



sx 



and 



X~ " X „ 

L 0 ~ 



1 

1 V v sx' Vo 



( 22 ) 



for the smaller and larger stationary points respectively. Eq. (21), giving 
the sensitivity to coupling for the two-way matched case (Case 2), is seen 
to be identical to the sensitivity for Case 1, the one-way coupling condition 
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given in Eq (11). Eq (22), shows that the second stability point is always 

very large, i.e. much larger than one and therefore lies outside the 

physically reliazable uncertainty space. 

The important thing to conclude from this analysis is that for small 

counter-C coefficients, the relative loss of knowledge is proportional 

3 

to the ratio of counter-C effort to the target's natural entropic birth 
rate. If the target's information decay is very slow under normal conditions, 

3 

then he will be very susceptible to counter-C activities. However, if the 
normal environmental relaxation times, T , are short, (his entropic birth 

rate is high), more effort will be required to increase the target's average 

\ 

uncertainty an equivalent proportion. 

V .2 Stability 

The preceding paragraph dealt strictly with stationary or equilibrium 
behavior and analyzed the sensitivity of the equilibrium condition to varia- 

3 

tion of system parameters, in particular the counter-C coefficients 

y & y . In this section, we wish to investigate the dynamic behavior 
yx xy 

of the system. We can describe rather thoroughly the system behavior of 

this non-linear system for small deviations from equilibrium using a tech- 

(4) 

nique from non-linear mechanics. 

The analysis begins by translating the equations to the stationary 
point by the change of variables 

x' = x-x-| (23) 

y' = y-y-i 

which leads to 

x' = t( ot 2 y l ' a l^ x ' " a 2 ^ 1_x l^' + a 2 x ' y 'J 
y' = [(b 2 x 1 - 3-, )y ' - 3 2 (.l-y 1 )x' + e 2 y'x'] 

(4) 

v 'See, for example; Minorsky, "Non-Linear Mechanics", J.W. Edwards, 

Ann Arbor, 1947. 



max 



max 



(24) 
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For small displacements from equilibrium, the product terms may be 
neglected so that motion near the stationary point is described by the 
solution of the coupled linear dynamical equations, 
x'= [a^ - ct-| ) x * - a 2 (- Xl )y'] Y max 



( 25 ) 



y<= [(B 2 x r 3-j )y 1 - B 2 0-y 1 )x'] X max 

Before undertaking a general study of these equations, let us consider 
conditions analogous to those of Case 1 described in part V.l above, i.e. 
one-way countering of Y on X such that S 2 =0. For simplicity, we shall 
further prescribe that X max = Y max , and that v $x + v ix = v $y + v. y . With 
these restrictions we are led to solutions 
y 1 = y' (o) exp [ - (V IY + V SY )t] 



x‘ = x' (0)exp [- (V 
where x'(0) and 
V I Y * V SY’ V IX’ V SX 



IX tv SX + V (1 ' X l )X max )t l- y 'rR^T l 1 - exp t- Y yx {1 ' X l ) Wl] 
y'(0) are the initial perturbations and we recall that 
are the non-normal i zed birth and death rates for X & Y. 



We see that Y, which is not being actively countered, has the same 

dynamic behavior as in Eq (3), simple exponential decay back to equilibrium 

-1 1 s 

with time constant (V + V. ) However, perturbations of X entropy 

sy i y 

are actually forced back to equilibrium more rapidly. Moreover, although 
perturbations in X do not affect Y, displacements of Y from equilibrium 
do cause variations in X, but always in the opposite sense . Thus, if Y 
has a temporary loss of knowledge, X will obtain a temporary increase, and 

I c 

vice versa (See Figure 4). Note that X maximum good fortune (or bad as 
the case may be) will be delayed from the time of Y maximum loss of 
knowledge. The time delay of the maximum is given by 




In 



v 1x + v sx 



V Ux l ,/(v i'c +v sx> 



(27) 
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“ here t r = (V IX + v sx) 



•1 



is the uncoupled relaxation time of the system. 

Eq (27) is plotted in Figure 5. We see that for very small couplings, the 
maximum effect is delayed about one system time constant, whereas for 
large couplings, the maximum effect occurs almost immediately. 

A more general analysis of Eq(s) (25) is made easier by considering the 
reversible linear transformation 

C = A x 1 + By 1 (28) 

n = Cx' + Dy' 

with A, 8 , 0,818 chosen such that 

k = 



% = 5 e 

^ ^0 



V 



n = s 2 n 



n = n 0 e 



s 2 t 



(29) 



and 



A B 
C D 



j* 0. Eq's (29) are known as the cannonical form, 



(4) 

Following Minorsky , one finds that the exponential coefficients. 



S-j 2 are given by roots of the characteristic equation, 

S 2 - p S + q = 0 
S 1 ,2 = P/2 + \ (p 2 -4q) % 
where in our case: 

p = {a 2 y 1 - ot 1 ) Y max + (8 2 x 1 - 3 1 )X max 

" ' “2 e 2 - V<fj - xp-o-xpn-y,)] Y max 

p 2 -4q - [<Vl- a l )Y max ' (6 2 x r 6 l )X max] 2 + 4a 2 6 2 (, - x l > (1 'X, > Wmax 



(30) 



(31) 



The nature of the behavior of the cannonical equations, and hence 
x‘ and y' through the inverse linear transformation, is determined by the 
location in p,q space. Figure 6 shows the type of dynamical behavior that 
obtains near equilibrium in various regions of p,q space. 
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It is also obvious that 



Since y- | and B-j > $ 2^1 , p < 0. 

2 

> 0 so that P _4< 1 > 0 and we are never in Regions II or III of 

Figure 6. We can also see that since ^1_ >0, and ^1_ >1, that 

<*2 &o 

q > 0 since CI 2 B 2 > 0. 

Therefore, we see that p and q are always found in Region IV of 
Figure 6 and the equations exhibit stable nodal point behavior near 
equilibrium. That is to say, in Region IV, S-j and S 2 are always real 
negative numbers, and therefore E, and n are simple damped exponentials. 

Our small perturbation solutions are then of the form (see Minorsky, 
P. 44) 



s t 

x ' = [ S 2 + Y max (c, l - a 2*l>] e ' + [ S 1 + Y n,ax (a r <*2*1 >] % e 



S 2 t 



(32) 



S,t Sot 



X iax ¥ 1 - Y l>5o e ' + n o e 



which are just linear combinations of the damped exponentials. Thus the 

- 1 -1 

system is returned to equilibrium with the time constants -S-| and -S 2 
determined according to Eq's (30). The constants and n Q satisfy the 
initial conditions 

x 1 (0) = (S 2 + a-| - a 2 y-|)C 0 + (S-j + a-j - a 2 y-|)n o 



y'(o) = B 2 (l-y 1 )(? 0 + n 0 ) 



(33) 



Consider now the nature of these solutions for conditions somewhat 
analogous to Case 2 of part V.l, the "matched case". In particular, again 

,et X max“ Y max' V Y xy’ (ct 2' 6 2> ' V*1 and v sx'V (This assure 
that a^= B-j). Under these "dynamically matched" system conditions we 

find that 



p = 2 (a 2 x 1 - a 1 )X max 
ls> U,)h , 2a 2 (l-x 1 )X max 



(34) 
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so that the characteristic roots of Eq (30) are simply 



5 1 " - ( a i" a 2^ X max = - ^ V IX + V SX^ 

5 2 = S 1 ' 2 Y yx (1 " x l )X max 

The dynamic behavior (Eq's (32)) becomes 



( 35 ) 



S,t S n t 

x' = ( i ) e 2 + (A ) e 1 



S 0 t 



Sit 



(36) 



(1) e 2 - (A) e 1 



where, 5* = x'(0) + y'(0) , A = x'(0) - y'(0) 



Suppose the system is initially disturbed so that x'(0) = -y'(0). 
Then Eq's (36) become 



x' = x'(0) e 



" (V IX + V SX^ 



y' = y'(0) e 



“( V IY + V SY^ 



(37) 



which are identical to a completely uncoupled system (See Eq (3)). However, 
suppose the system receives perturbations of the same signs, i.e. 
x'(0) - y'(0). Then 



x' 

y' 



x 1 (0) 

y'( o) 



e -< v ix + v sx + V'VW 1 

-< V IY + V SY * V'l’W' 
e 



(38) 



and the systems are driven back toward equilibrium even faster than when 

3 

no coupling exists by virtue of the counter-C activities. Finally, 
suppose x'(0) = x'(0) and y'(0) = 0 . Then 



x' = 



x'(0) 

~ T~ 



r 



--x'(O) 

2 



e '^ V IX + V IY )1: ^ + e " 2 Y yx ( 1 " x l^ X max t ] 
e -(Vjx + v iY)t ^ _ e -2 Y xy( 1 -y 1 )v max t] 



( 39 ) 
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and Y will experience a change of entropy in the opposite sense of X . 



/ 



Again the maximum effect will be delayed as in the one way coupling case 

3 

according to Eq (27), but with the effect of the counter-C coefficient, 
y , doubled. 

Ay 

To summarize, the equilibrium points are always stable, and non- 

3 

oscillatory for the counter-C model investigated in this section. Small 
perturbations of the system from equilibrium are powerfully forced back 
to equilibrium with time constants equal to or less than those of the 

3 

individual systems. However, counter-C efforts do cause a perturbation 
in one system to appear as a perturbation in the other system, but with 
some time delay. 

V .3 Dynamic Behavior Far From Equilibrium 

It is difficult to find general results for the behavior of the system 
far from equilibrium. However, a modest numerical investigation of Eq's (24) 
has revealed the following: 

1) For initial perturbations of about 0.1 (i.e. one tenth the entire 
entropy range), the behavior seems well described by Eq's (32), the 
small perturbation result. 

2) For large initial displacements, e.g. 0.9, the behavior is still 
stable nodal point in character and the entropies return to equilibrium 
smoothly. However, the "modeled entropy" for one system may temporarily 
go negative, which is physically impossible, although neither entropy 
ever appears to exceed its maximum value. 

3) The rate of return to equilibrium seems proportional to the counter- 

3 3 

C coefficients. Strong counter-C activities appear to strengthen the 
forces returning the system to equilibrium (Also see Eq's (39)). 

We have shown some typical dynamic trajectories in Figure 7 for 2 
equilibrium points and 3 initial conditions. Note the asymmetry of the 
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trajectories returning to (.5,. 5) from (1,0) and (1,1). Also the one 
returning from (1,1), total uncertainty is restored to equilibrium much 
more quickly than the one coming from total uncertainty of X and perfect 
knowledge of Y, (same signs versus opposite signs). 

Two of the trajectories become negative although beginning and ending 
inside the physically realizable region. What will happen if a non-negativity 
constraint is imposed is not known, but could rather easily be investigated 
numerically. What does seem to be true is that even starting on a boundary, 
the entropies do not blow up or trend toward the other equilibrium point in 
the first quadrant but outside the physical region. One should be cautioned 
that this conclusion is based on a limited numerical investigation of the 
non-linear dynamical system equations and further investigation could indeed 
prove otherwise. 

VI . Discussion 

3 

The model we have proposed for C information dynamics based on entropy 

3 

includes a non-linear term to account for counter-C activities of two sides 
engaged in an information war. Before we discuss the meaning and implication 
of our analytical results, let us briefly recount some other applications of 
this class of equations. 

VI .1 Population Dynamics 

An interesting example is provided by the Lotka-Vol terra model of 
population dynamics, originally devised to explain temporal oscillations 
in the occurrence of fish in the Adriatic Sea. If we let the prey fish 
be x and the predator fish be y, then the model is described by, 

m 

x = a-jX - (*2 xy 

(40) 

y = -3-j y + B 2 xy 
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where and are the prey hirth and predator death rates respectively 

and 0.2 and $2 describe the prey losses and predator gains due to predator 
feeding. These equations have been studied extensively (along with more 
elegant versions) for describing population dynamics of conflicting species. ^ 
They have non-zero equilibrium points about which under some conditions they 
exhibit stable focal point behavior, unlike our counter-C^ equations. Note 
that the non-linear terms are of opposite sign. As we shall see shortly, 
this is a requirement for focal point behavior. 



VI .2 Models of Combat 

Another example is provided by the Lanchester equations for armed 
combat. ^ Here we take x and y to be the sizes of the forces of two 
sides. Then 



x = -a-jy - a 2 xy 
y = -$ 1 x -$ 2 xy 



( 41 ) 



include both terms most commonly considered. When a 2 = $2 = Eq's (41) 
are of the form Lanchester termed "aimed fire" or "modern warfare". When 
a-j = $i * 0, Eq's (41) are of the form he called area fire or "ancient 
warfare". Although a non-zero equilibrium point exists in the general case 

of Eq's (41), we see that no non-zero equilibrium is possible for either of 

Ax 

the special cases. However, an "exchange ratio", jy , is obtained by 
dividing the two equations for these two cases. The solutions show the 
trajectories that x and y must follow as the forces are depleted. In 
particular, for "aimed" fire, the exchange ratio is 



May, R.M., "Stability and Complexity in Model Ecosystems", 
Princeton Univ. Press, 1973... 

^Taylor, J.6., "A Tutorial on Lanchester-Type Models of Warfare", 
Proc of 35th MORS Symposium, July 1975. 
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X 

y 



( 42 ) 



dx 

dy 



a. 



3, 



e,tx? - x 2 ) = a,(y* - y 2 ) 



and for "area" fire 
dx _ ^2 

dy " 3 2 

3 2 (x. - x) = a 2 (y. - y) 



(43) 



These results have been interpreted to show the advantage in concentra- 
ting forces to minimize losses when "aimed" fire prevails whereas concentra- 
ting forces in "area" fire is of no particular advantage. To illustrate, 
let a 2 = 6 2 = a-j = 6-j =1 and let y^ = 0. Then for "area" fire 

x i - x f < s 'c 

= 1 regardless of x.. so X and Y losses are matched. 

■M 

However, for "aimed" fire 




which for — ^ - « 1 is approximated by 
x i 



x i ~ x f 1_ (45) 

y i * 2 X 1 

For example, if x^ = y^ then of course the losses are matched (Eq. (44)), 
however, if X initially concentrates all his forces against Y, so, say 

i c * c 

x.j = 4y.. , then X losses are only 1/8 of Y losses. 

One of the characteristics of Lanchester's specialized equations are 
that their dynamical solutions, like our C equations, are sums of real 
exponentials. However, since they have no non-zero stable point, their 
equilibrium behavior is not of great interest. Also, the solutions can go 
negative, which is not physically consistent with reality. 
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3 

VI. 3 Alternative C Information Models 



Consider the pair of equations 
x = a - a,x + f(x,y) 

(46) 

y = 6 0 ' e i y + s(y.x) 



Let us list some of the other options for f(x,y) and g(y,x), along 
with the one analyzed in Part V, and their interpretations. 



Alternative 


f(x,y) 


g(y.x) 


Entropy Birth (Death) 
Rate 


a.) 


a 2 (l-y) 


3 2 (1-x) 


3 

Counter-C efforts create 
entropy in proportion to 
own knowledge. 


b.) 


a 2 (l-x)(l-y) 


b 2 ( i -y) ( i -x) 


3 

Counter-C efforts create 
entropy in proportion to 
product of both sides 
knowledge. 


c.) 


a 2 x(l-y) 


6 2 y(l-x) 


3 

Counter-C efforts create 
entropy in proportion to 
self knowledge and opposi- 
tion's ignorance. 


d.) 


-a 2 y(l-x) 


-0 2 x( 1 -y) 


Intelligence efforts destroy 
entropy in proportion to 
self knowledge and opposi- 
tion's ignorance. 



2 

Alternatives a.), b.) and c.) are identical to the counter-C coupling 
options a.), b.) and c.) of Part IV. Alternative b.) is the case analyzed 
extensively in Part V. Alternative d.) listed here is a new concept. It 
proposes that instead of counter-C activities to increase the opponent's 
entropy, the players focus on special intelligence activities to reduce their 
own uncertainty and that the entropy death rate by such means will be pro- 
portional to one's own knowledge and the opposition's ignorance. 

The implications of these four alternatives for dynamic stability 
are listed below: 
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A1 ternative 


Re [Si & S 2 J 


Im[S 1 & S 2 J 


a.) 


Both < 0 


0 




(stable) 


(no focal points) 


b.) 


Both < 0 


0 


(stable) 


(no focal points) 


c.) 


Either or both 
may be > 0 
(may be unstable) 


0 

(no focal points) 


d.) 


Either or both 


0 




may be > 0 
(may be unstable) 


(no focal points) 



We see that models c.) & d.) may lead to unstable behavior near equilibrium 

points that would cause one or both entropies to blow-up. However, no model 

produces focal point behavior near equilibirum. This is because in all 

models, a 2 S 2 > 0, which forces the imaginary parts of the roots of the 

characteristic equation to be zero. (This is not so with the Lotka-Vol terra 

equations for population dynamics.) 

An interesting avenue for further research will be to investigate 

combinations of counter-C and intelligence models both for stable and 

unstable as well as focal point behavior and for equilibrium sensitivity 

3 

to intelligence vice counter-C efforts. 
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VII . Summary 



3 

In this paper we have proposed a model for C information dynamics 

3 

incorporating the effects of counter-C activities. The model is based 
on the inevitable growth of uncertainty inherent in military situations 
and the concept of information sensors and sources acting as constraints 
that maintain uncertainty below its worst possible value. System entropy 
dynamics of two opposing sides are characterized by natural "birth 11 and 

3 

"death" rates of entropy. Counter-C activities are introduced as additional 
growth terms that depend in some way on the entropy of one or both players. 

The dependence analyzed extensively in this paper in Part V models 

3 

counter-C effectiveness as being in direct proportion to the product of 
the two system's knowledges, (X m , -X) (Y m= -Y), where X and Y are the 
two entropies and X_ 3V and Y m=v the largest amount of uncertainty (or 
information) possible in the two systems respectively. It is shown that 
for this kind of coupling, the relative shift of system equilibrium is 
directly proportional to the ratio of the coupling coefficient to the 
system's natural entropy birth rate. Furthermore , it is shown that small 
perturbations from equilibrium are restored to equilibrium by the system 
forces, i.e. the system is ultra-stable, but that 1) perturbations of X 
and Y with the same sign are restroed much more rapidly than perturbations 
of opposite sign, and 2) a perturbation in one system induces a delayed 
perturbation in the other system of the opposite sign. Thus, if X becomes 
fortuitously more knowledgeable by chance, Y will in turn some time later 
become more ignorant and vice versa. It is also shown the system's dynamical 
trajectories near equilibrium are described by sums of exponentials with real 
coefficients. Such equilibrium points are called "nodal points". (In 
contrast, "focal point" systems have exponentials with complex coefficients 



29 



in which the trajectories spiral into the equilibrium point. See Figure 6.) 

A modest investigation of the dynamics far from equilibrium show the 
system always returns to equilibrium; but it is possible for the "modeled 
entropy" of one or both systems to temporarily become negative, a physically 
impossible condition for the "real entropy". The importance of a non- 
negativity constraint on system behavior far from equilibrium requires 
further investigation. 

Finally, after a brief review of two other well known applications of 
coupled non-linear state equations, the Lanchester combat equations and the 

O 

prey-predator equations, several alternative counter-C models and an 

"intelligence" model are proposed. It is shorn that these models too exhibit 

nodal point behavior, unlike for example, the predator-prey model, which may 

exhibit focal point behavior. However, the intelligence model and a 
3 

counter-C model where the entropy rates depend on the product of one's own 
knowledge and the opposition's ignorance, are not necessarily stable. That 
is, small deviations from equilibrium may cause one, or the other, or both 
entropies to diverge. 

It is clear that further theoretical investigation of entropy models 

is needed, along with some modest simulations of actual systems, to determine 

3 

the utility of this approach for the evaluation of specific C I alternatives. 

NPS 

Monterey, CA 
December 1980 
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